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Abstract 

LISA  ( Laser  Interferometer  Space  Antenna)  is  a  proposed  space  mission  which  will  use  coherent 
laser  beams  exchanged  between  three  remote  spacecraft,  to  detect  and  study  low-frequency  cosmic 
gravitational  radiation .fl-1  The  multiple  Doppler  readouts  available  with  LISA ,  which  incorporate 
frequency  standards  for  measuring  phase  differences  between  the  received  and  transmitted  laser 
beams,  permit  simultaneous  formation  of  several  observables  J234!  All  are  independent  of  lasers  and 
frequency  standard  phase  fluctuations,  but  have  different  couplings  to  gravitational  waves  and  to 
the  various  LISA  instrumental  noises .  Comparison  of  the  conventional  Michelson  interferometer 
observable  with  the  fully -symmetric  Sagnac  data-type  allows  unambiguous  discrimination  between  a 
gravitational  wave  background  and  instrumental  noise .  The  method  presented  here  can  be  used  to 
detect  a  confusion-limited  gravitational  wave  background . 


INTRODUCTION 

The  Laser  Interferometer  Space  Antenna  (LISA)  is  a  space  mission,  jointly  proposed  by 
NASA  and  ESA,  aimed  to  detect  and  study  gravitational  radiation  in  the  millihertz  fre¬ 
quency  band.  With  its  three  spacecraft,  each  carrying  lasers,  beam  splitters,  photodetectors 
and  drag-free  proof  masses  on  each  of  their  two  optical  benches,  LISA  will  have  the  ca¬ 
pability  of  measuring  six  time  series  of  Doppler  shifts  of  the  one-way  laser  beams  between 
spacecraft  pairs,  and  six  shifts  between  adjacent  optical  benches  on  each  spacecraft.  By 
linearly  combining,  with  suitable  time  delays,  these  twelve  data  sets,  it  will  be  possible 
to  cancel  the  otherwise  overwhelming  phase  noise  of  the  lasers  (Av/v  ~  10-13)  to  a  level 
h  ~  A v/c  ~  10-23.  This  level  is  set  by  the  buffeting  of  the  drag-free  proof  masses  inside 
each  optical  bench,  and  by  the  shot  noise  at  the  photodetectors^. 

LISA  is  expected  to  detect  monochromatic  radiation  emitted  by  galactic  binary  systems. 
Particularly  at  low  Fourier  frequencies  (say  0.1  —  8  mHz),  however,  there  will  be  many 
galactic  binaries  radiating  within  each  Fourier  resolution  bin^h  These  latter  signals  will  not 
be  detectable  individually,  forming  a  continuum  which  could  be  confused  with  instrumental 
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noise.  The  level  of  this  stochastic  background  is  uncertain,  but  could  be  in  the  range  10-20  — 
10~23.  Since  these  galactic  binary  populations  are  virtually  guaranteed,  the  detection  of  their 
signals  could  be  the  first  direct  detection  of  gravitational  waves. 

For  this  measurement  it  is  very  desirable  that  competing  proof-mass  and  other  instrumen¬ 
tal  noises  be  both  characterized  and  calibrated  before  flight,  and  measured  in  the  actual  flight 
configuration  while  data  are  being  taken.  In  contrast  to  Earth-based,  equal-arm  interferom¬ 
eter  detectors  of  gravitational  radiation,  LISA  will  have  multiple  readouts,  and  the  Doppler 
data  they  generate  can  be  combined  differently  to  give  measurements  not  only  insensitive 
to  laser  phase  fluctuations  and  optical  bench  motions,  but  also  with  different  sensitivities  to 
gravitational  waves  and  to  the  remaining  system  noisef3,45. 

In  this  .article  we  discuss  two  laser-and-optical-bench-noise-free  combinations  of  the  LISA 
readouts,  previously  denoted  (  (Sagnac)  and  X  (Michelson),  that  have  very  different  re¬ 
sponses  to  the  gravitational  wave  background  but  comparable  responses  to  instrumental 
noise  sources^4’5!. 


THE  SAGNAC  AND  MICHELSON  INTERFEROMETERS 


The  six  Doppler  beams  exchanged  between  the  LISA  spacecraft  imply  the  six  Doppler  read¬ 
outs  yij  ( i,j  =  1,2,3)  recorded  when  each  transmitted  beam  is  mixed  with  the  laser  light 
at  the  receiving  optical  bench.  Delay  times  for  light  travel  between  the  spacecraft  must  be 
carefully  accounted  for  when  combining  these  data.  Six  further  data  streams,  denoted  Zij 
(i,j  =  1,2,3),  are  generated  internally  to  monitor  both  lack  of  rigidity  and  laser  synchro¬ 
nization  between  the  independent  optical  benches  at  each  spacecraft.  The  combination  ( 
uses  all  the  Doppler  data  symmetrically^4’5] 

C  =  2/32,2  —  2/23,3  +  2/13,3  —  2/31,1  +  2/21,1  —  2/12,2 

1 

+  2  (“213,21  +  ^23,12  —  321,23  +  ^31,23  ~  ^32,13  +  ^12,13) 

+  2  (~Z32>2  ■*“  Zl2>2  ~  Zl3>3  Z23>3  ~  321,1  +  331, l)  ■  (1) 

The  comma  notation  indicates  time-delays  along  the  arms  of  the  3-spacecraft  configuration 

2/32,2  =  2/32 {t  —  L2)  ,  (2) 

and  so  forth  (units  in  which  c  =  1). 

The  transfer  functions  of  (  to  instrumental  noises  and  to  gravitational  waves  were  cal¬ 
culated  in  references  [3,  4].  The  resulting  instrumental  noise  power  spectrum  for  £  is  shown 
in  Figure  1.  Also  shown  there  is  the  computed  power  spectrum  of X,  averaged  over  the  sky 
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and  elliptical  polarization  states,  that  would  result  from  a  stochastic  background  originated 
by  an  ensemble  of  galactic  binary  systems^ . 

The  laser-and  optical-bench-noise-free  combination,  X,  only  requires  four  data  streams. 
This  combination  is  equivalent  to  an  (unequal  arm)  Michelson  interferometer.  Its  expression 
is  equal  to^4*5! 


X  =  2/32,322  —  2/23,233  +  2/31,22  ~  2/21,33  +  2/23,2  —  2/32,3  +  2/21  ~  2/31 

+  2^_ Z21>2233  ^21,33  +  -221,22  —  z2l) 

1 

+  “(+-231,2233  —  -231,33  —  -231,22  +  -23l)  •  (3) 

The  expected  instrumental  noise  power  spectrum  in  X  is  shown  in  Figure  1.  Also  shown 
is  the  anticipated  galactic  binary  confusion  spectrum)1!,  which  would  be  observed  in  X. 
Comparison  of  X  and  (  allows  the  background  to  be  discriminated  from  instrumental  noise. 


DETECTING  THE  GALACTIC  STOCHASTIC  BACKGROUND 


The  flight  configuration  of  the  three  spacecraft  forming  LISA  will  be  essentially  equilateral, 
with  Li  =  L2  =  L3  =  L  =  16.67  sec.  In  the  frequency  band  of  interest  Jfi.l  —  8  mHz),  the 
expressions  for  the  Fourier  transforms  of  the  gravitational  wave  signals  Xgw(f),  ( 9W(f)  and 
the  power  spectral  densities  of  the  system  noises  in  X  and  (  ?  S xn°i3e  (/),  sc  noise  (/),  can  be 
Taylor-expanded  in  the  dimensionless  quantity  fL.  The  first  non-zero  terms  are  equal  to 


X9W(f ) 
Cgw(f) 


S Xn°i' se  (/) 


S^noise  (/) 


2  (2m fL)2  [n3  •  h(/)  •  h3  -  n2  ■  h  (/)  •  n2  , 

^  (2mfLf  [(^  •  hi)(hi  ■  h (/)  •  hi)  +  (k  ■  h2)(h2  ■  h (/)  •  h2) 
+  (k  -n3)(n3  •  h(/)  -n3)]  , 


S Xvroof-mass  (/)  +  Sxopticalpath(f) 


16  [S,  (/)  +  Sr  (/)  +  S3(f)  +  S2.  (/)]  (2tt !Lf 
+4[<S32(/)  +  S2s(f )  +  53i(/)  +  S2i(f)\  (27 r/L)2 
[Si(f)  +  S2(f)  +  S3(f)  +  Sr  (f)  +  s2.  (/)  +  s3.  (/)]  (27 r/L)2 
+  [<^32 (/)  +  S23(f)  +  S3i(f)  +  S2i(f)  +  Si3(f)  +  S12(f)}  , 


(4) 

(5) 

(6) 
(7) 


where  we  have  denoted  by  SXvroofmass(f ),  and  SXopticabath(f )  the  aggregate  contributions  to 
the  power  spectrum  of  the  noise  in  the  response  X  from  the  proof  mass  and  optical  path 
noises  respectively.  The  expressions  in  square  brackets  in  Equations  (4,  5)  incorporate  LISA 
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antenna  responses f3,4^,  and  are  of  the  same  order  of  magnitude.  The  proof  mass  Doppler 
noise  spectra  5j(/),  Si*(f )  ( i  =  1,2,3)  will  be  designed  to  a  nominal  power  spectral  levelM 
S°(f)  =  2.5  x  10~48[f/lHz]~2  Hz-1,  while  the  optical  path  noise  spectra  Sij(f )  ,(i,j  = 
1,2,3  ,i  t -  j),  which  include  shot  noises  at  the  photo  detectors  and  beam  pointing  noise W, 
are  expected  to  be  equal  to  a  nominal  spectrum  S1(f)  =  1.8  x  10~87[f /1Hz]2.  Both  these 
noise  sources  will  be  estimated  before  launch,  but  could  be  larger  when  the  in-orbit  data  will 
be  taken. 

First  consider  the  responses  to  the  gravitational  wave  signal,  given  in  Equations  (4,  5). 
At  /  =  10-3  Hz,  for  instance,  (where  27 r/L  ~  10-1)  the  absolute  value  of  the  coefficient  in 
front  of  the  squared-bracket  in  the  £  response  (Eq.  5)  is  about  three  orders  of  magnitudes 
smaller  than  the  corresponding  coefficient  given  in  the  expression  for  X  (Eq.  4).  The  power 
spectral  densities  of  the  noises  due  to  the  proof  masses  and  the  optical-path  noise  (Eqs.  6,  7) 
will  only  differ  at  most  by  an  order  of  magnitude.  We  conclude  that  in  this  lower  frequency 
range  the  LISA  Sagnac  response,  £,  can  be  used  as  a  gravitational  wave  shield.  In  what 
follows  we  will  ignore  the  gravitational  wave  background  contribution  to  £. 

To  take  quantitative  advantage  of  this  property  of  £,  consider  the  observed  power  spectral 
densities  of  X  and  £ 


Sf(f) 

S?(f) 


SxAf)  +  s  Xproofmasa  (f')  td  ^optica  Ipatk  (,/') 


1 

16 


7  V  proo  fmas 


s(/)  + 


SXopt 


icalpath 


(/) 


WP2 

+  [S2(/)  +  S3-(/)]  (2<r /L)\ 


+  [Sl3(/)  +  S12(/)] 


(8) 

(9) 


where  in  Equation  (9)  we  have  written  the  power  spectra  of  the  noises  in  £  in  terms  of  the 
power  spectra  of  the  noises  in  X  and  of  some  remaining  terms  that  are  not  present  in  X, 
to  emphasize  commonality  of  some  noise  sources.  We  suppose  that  the  noise  contributed 
by  any  one  of  the  proof  masses  and  optical-path  noise  sources  will  be  greater  than  or  equal 
to  the  design  values,  S°(f)  and  Sl(f)  respectively.  From  Equation  (9),  if  the  magnitude  of 
the  measured  power  spectral  density  of  the  response  £  is  at  its  anticipated  level  S°bs(f  )  — 
6  5°(/)(27t/ L)2  4-  6  S1(f),  then  the  level  of  the  power  spectral  density  of  the  noise  entering 
into  X  is  known.  The  spectrum 

SxAf)  =  S$s(f)  -  64  S° (f)(2ir f L)2  -  16  51(/)(2tt/L )2  ,  (10) 


should  then  be  attributed  to  a  galactic  binary  background  of  gravitational  radiation.  In  any 
event,  the  RHS  of  Equation  (10)  is  an  upper  bound  to  Sxs »■ 

On  the  other  hand,  if  the  measured  spectrum  of  £  is  above  its  anticipated  design  level, 
consider  the  following  combination  of  the  measured  spectra 
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(11) 


St  (/)  -  16  Sf(f)  = 


5xr~'  -  16 '.S'jf/j  -T  Sy  ( f  ]\(2x  f  L'f 

-16[S13(/)  +  Sb(/)] 

—16  [532(/)  +  S'xJ.  f )  +  SVi  if)  +  S2  ! /)]  x 
x  [l  -  (t r/L)21 


The  coefficient  of  S£bs  has  been  chosen  so  that  the  noise  terms  on  the  right  hand  side  are  all 
now  negative-definite  and  can  thus  be  bounded  from  above  by  their  design,  or  nominal,  values 
S°(f)  and  S1(f)  respectively.  The  result  is  a  lower  bound  for  observational  discrimination 
of  the  gravitational  wave  background  spectrum 

Sx,-U)  >  ££“(/)- 16  Sc"l,(/)+ 32  (27t/L)2S°(/) 

+16  [6  -  (2 xfL)2}S\f).  (12) 

Equations  similar  to  (11)  and  (12)  can  be  written  for  the  other  two  interferometer  combina¬ 
tions,  Y  and  ZM.  In  those  equations,  there  will  be  different  mixes  of  canceled  and  bounded 
noise  sources,  resulting,  in  general,  in  different  gravitational  wave  spectrum  lower  bounds. 


CONCLUSIONS 

The  response  of  the  Sagnac  interferometer  to  a  gravitational  wave  signal  is  several  orders 
of  magnitudes  smaller  than  that  of  the  Michelson  interferometer.  In  the  frequency  band  of 
interests  (0.1  —  8)  mHz,  however,  the  Sagnac  response  to  the  noise  sources  is  of  the  same 
order  of  magnitude  as  that  of  the  Michelson  interferometer.  As  a  consequence  of  these  facts 
we  have  shown  that  it  is  possible  to  estimate  the  magnitude  of  the  noise  sources  affecting 
the  Michelson  interferometer  response  in  the  low-frequency  region  of  the  accessible  band  by 
using  the  Sagnac  interferometer.  This  in  turn  allows  us  to  discriminate  a  gravitational  wave 
background  of  galactic  origin  from  instrumental  noise  affecting  the  Michelson  interferometer 
response. 
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Questions  and  Answers 


ROBERT  NELSON  (Satellite  Engineering  Research  Corporation):  It  might  be  appropriate  to 
note  that  Dr.  Joseph  Webber  of  the  University  of  Maryland  passed  away  this  past  September. 
Dr.  Webber  founded  the  whole  subject  of  experimental  gravitational  wave  detection  at  the 
University  of  Maryland.  In  addition,  in  the  late  1940s,  he  published  one  of  the  first  papers  on 
the  notion  of  a  population  inversion  concept  that  was  the  foundation  for  the  invention  of  the 
maser. 

Over  the  last  decade  or  so,  Dr.  Webber  did  analyses  that  indicate  some  evidence  of  gravitational 
wave  detection  by  VAR  detectors  at  the  University  of  Maryland  and  in  Rome  which  were 
correlated  to  observations  of  neutrinos  that  were  associated  with  the  supernova,  I  believe  it 
was  1985.  One  of  the  things  that  you  might  consider  in  the  future  is,  in  light  of  that  type  of 
analysis,  is  to  correlate  your  own  gravitational  wave  measurements  with  other  sources  such  as 
neutrinos  which,  according  to  theory,  are  produced  in  gravitational  collapse  events. 

MASSINO  TINTO:  Yes,  I  would  like  to  add  and  I  should  have  said  it  earlier,  this  particular 
instrument  will  work  in  the  millihertz  frequency  band.  So  we  are  expecting  to  observe  sources 
that  will  not  be  in  the  kilohertz  frequency  band  where  less  detectors  were  operating.  So 
supernova  light  is  frozen  and  will  not  very  likely  be  in  this  span  of  observations.  Usually 
supernova  explosions  are  accompanied  by  neutrino  emissions. 

The  reason  for  this  analysis,  in  a  sense,  is  really  trying  to  assess  where  the  sensitivity  is  without 
relying  on  other  instruments.  In  a  sense,  it  would  corroborate  your  observation.  So  once  you 
know  what  the  sensitivity  curve  is,  anything  that  is  above  it,  you  know,  sort  of  classifies  as  other 
sources. 

DEMETRIOS  MATSAKIS  (USNO):  Is  there  any  timing  requirements  that  you  could  talk  about 
with  this  system? 

TINTO:  Yes,  and  I  am  getting  to  that.  You  see,  when  you  actually  phase  for  the  phase  differences 
from  these  Doppler  measurements,  you  have  to  take  into  account  that  the  spacecraft  are  moving 
relative  to  each  other  and  is  not  stationary.  So  at  10~14  hertz,  which  is  the  frequency  of  the 
laser,  the  relative  speed  of  10  meters  per  second  introduces  big  nodes  of  several  megahertz. 
You  have  to  track  that  bit  node  in  your  phase  measurement.  To  do  that,  we  rely  on  USOs,  so 
we  have  a  timing  system  aboard  this  spacecraft  that  allows  us  to  measure  the  phase  and  then, 
on  the  other  hand,  introduce  noise.  It  is  so  noisy,  and  we  want  to  remove  that  noise.  There 
are  techniques  to  eventually  remove  this  noise  measurement  introduced  by  the  USOs  which 
are  incorporated  into  this  interferometric  technique  that  we  have.  So  the  system  is  certainly 
has  on-board  timing  systems  and  timing  requirements.  In  fact,  we  rely  on  state-of-art  USOs  in 
order  to  do  these  measurements.  We  need  parts  in  1015  or  better  USO  performances. 

But  I  didn’t  want  to  present  all  the  details  of  the  timing  system  on  board  the  spacecraft. 

MATSAKIS:  It  is  covered  in  your  published  paper  though. 

TINTO:  Yes,  this  is  published  and  will  be  in  the  Proceedings  anyway. 

THOMAS  CLARK  (NASA  Goddard  Space  Flight  Center):  I  had  looked  at  a  similar  type 
interferometer  at  one  time,  and  maybe  I  don’t  understand  the  measurement.  If  you  are  talking 
about  the  laser  measuring  carrier  phase  of  the  laser  — if  I  did  the  calculations  you  are  talking 
about,  travel  times  of  about  16  seconds  — that  means  that  you  have  to  have  an  intrinsic  oscillator 
stability  at  the  laser  frequency  to  something  better  than  one  sixteenth  of  a  hertz.  Probably 
more  like  10  millihertz  at  1014  hertz.  So  that  says  you  are  expecting  a  laser  to  have  an  intrinsic 
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frequency  stability  over  the  Allan  variance  at  the  travel  time,  16  seconds,  of  something  like  a 
part  in  1016. 

TINTO:  No,  no,  we  don’t  need  that.  Even  parts  in  1013  you  have  to  know  accuracy  of  the 
uplink.  The  separation  of  the  spacecraft  is  only  50  meters  or  so.  If  you  know  the  uplink,  then 
you  can  combine  to  six  measurements  in  such  a  way  as  to  remove  the  fluctuations  of  the  laser 
to  the  level  that  you  want  to. 

CLARK:  The  other  question  I  was  going  to  ask  has  to  do  with  the  background  noise,  which 
is  a  different  one.  I  will  have  to  think  about  your  comment  on  that  one.  When  we  had 
thought  about  this  in  terms  of  a  slightly  different  program,  rather  than  the  gravitational  wave 
background  noise,  the  gravity  field  background  noise  looked  to  me  to  be  a  serious  problem. 
For  instance,  the  gravity  field  changes  at  this  interferometer  location  due  to  Venus  orbiting 
around  Mars  are  many  times  the  effect  of  the  signal  that  you  are  seeing.  The  gravity  noise 
due  to  essentially  all  of  the  undetected  and  unmodeled  asteroids  looked  to  me  like  they  would 
give  so  much  position  noise  across  the  interferometer  that  gravity  field  could  not  be  detectable. 
Have  you  thought  in  terms  of  just  the  solar  system  gravity  noise? 

TINTO:  I  think  that  has  been  included  in  the  Phase  A  report,  which  you  will  find  at  that  Web 
site  that  I  gave.  I  believe  the  conclusion  was  that  that  was  not  a  problem.  Also,  you  have 
to  take  into  account  the  time  scales  that  you  are  talking  about.  LISA  will  be  1,000  to  10,000 
seconds  or  so  up  to  1  second  integration  time.  Most  of  the  facts  we’re  concerned  about  is 
the  time  scale  completely  outside  the  band.  It  will  not  be  affecting  the  performance  of  the 
instrument.  Other  things  have  been  analyzed  and  found  not  to  be  significant.  To  get  the 
details,  I  can  point  out  the  references. 

But  I  think  the  key  point  here  for  LISA  is  that  you  can  actually  synthesize  an  interferometer, 
which  means  you  can  remove  the  fluctuations  of  the  laser.  It  requires  a  knowledge  of  the 
uplink.  So  going  to  your  first  question,  if  you  know  the  uplink  with  an  aperture  of  30  meters 
or  something,  then  you  can  combine  these  six  Doppler  measurements  in  such  a  way  as  to 
remove  those  fluctuations.  So  you  need  the  laser  at  part  in  the  1013  or  so.  Even  worse,  if  your 
knowledge  of  the  uplink  is  better. 
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